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Abstract 

We present a precise measurement of the standard model CP-violation parameter sin 201, the 
direct CP violation parameter |A|, the lifetimes of charged and neutral B mesons and their ratio, 
and the B^-B^ mixing parameter Am^ based on a sample of 152 x 10^ BB pairs collected at the 
T(4S') resonance with the Belle detector at the KEKB asymmetric-energy e^e~ collider. One of 
two B mesons is fully reconstructed in a CP-eigenstatc or a flavor-eigcnstatc decay channel. The 
flavor of the accompanying B meson is identified from its decay products. From the distributions of 
the time interval between the two B meson decay points, we obtain sin2</)i = 0.728 ± 0.056(stat) ± 
0.023(syst), |A| = 1.007 ± 0.041 (stat) ± 0.033(syst), tbo = [1.534 ± 0.008(stat) ± O.OlO(syst)] ps, 
tb+ = [1.635 ± O.Oll(stat) ± O.Oll(syst)] ps, tb+/tbo = 1.066 ± 0.008(stat) ± 0.008(syst) and 
Am^ = [0.511 ± 0.005(stat) ± 0.006(syst)] ps~"^. The results for sin20i and |A| are consistent with 
the standard model expectations. The significance of the observed deviation from unity in the 
lifetime ratio exceeds five standard deviations. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 



In the standard model (SM), CP violation arises from an irreducible phase in the weak 
interaction quark-mixing matrix [Cabibbo-Kobayashi-Maskawa (CKM) matrix] [1]. In par- 
ticular, the SM predicts a CP-violating asymmetry in the time-dependent rates for and 

decays to a common CP cigenstate fcp-, where the transition is dominated by the h ccs 
process, with negligible corrections from strong interactions [2]: 

_ r[B°(t) ^ fcp] - r[B°(f) ^ /cp| „ 
-^f*) = rp«(i)^/.p] + r[BO(t)^/cH = 

where T[B^{t), B^{t) fcp] is the rate for 5° or B° decay to fcp at a proper time t after 
production, ^/ is the CP eigenvalue of fcp, Am^ is the mass difference between the two B^ 
mass cigcnstatcs, and (pi is one of the three interior angles of the CKM unitarity triangle, 
defined as 0i = vr — SiYg(y^lVtd/V*iycd)- Non-zero values for sin2(^i have been reported by 
the Belle and BaBar collaborations [3-5]. Belle's latest published measurement of sin20i 
is based on a 78 fb~^ data sample (data set I) containing 85 x 10^ BB pairs produced at 
the T(4S') resonance. In this paper, we report an improved measurement incorporating an 
additional 62 fb~^ (data set II) for a total of 140 fb~^ (152 x 10^ BB pairs). A precise 
knowledge of sin 20i is essential for testing the Kobayashi-Maskawa model of CP violation. 

The sin 201 measurement requires a determination of a proper-time resolution function 
and of the wrong-tag fractions using a large sample of exclusively reconstructed flavor- 
eigenstate decays. We perform a precise measurement of the mixing parameter Am^ and 
of the neutral (charged) B meson lifetime r^o (r^+) as a byproduct of this procedure. Our 
previous results are based on a 29.1 fb^^ data sample [6-8]; thus our new measurements 
with a 140 fb~^ data sample provide significant improvements. 

Changes exist in the analysis with respect to OTir earlier results. We apply a new 
proper-time resolution function that reduces systematic uncertainties. We introduce b- 
fiavor- dependent wrong-tag fractions to accommodate possible differences between B^ and 
S° decays. We also adopt a multi-parameter fit to the fiavor-eigenstate samples to obtain 
Am^, Tbo, Tb+, the resolution parameters and wrong-tag fractions simultaneously. There 
are other improvements in the estimation of background components that are made possible 
by the increased statistics. 

The data were collected with the Belle detector [9] at the KEKB asymmetric-energy e^e" 
collider [10], which coUides 8.0 GeV e~ on 3.5 GeV e+ at a small (±11 mrad) crossing angle. 
We use events where one of the B mesons decays to fcp at time tcp, and the other decays 
to a self-tagging state /tag, which distinguishes 5° from B^, at time ttag- The CP violation 
manifests itself as an asymmetry A{At), where At is the proper time interval between the 
two decays: At = tcp — ^tag- We also use events in which fcp is replaced by a fiavor 
eigenstate /av; the decay chain in this case is T(45') B^B^ — >■ /flv/tag- The time evolution 
is described as e^l^*l/^B" /(4r5()){l ± cos(Am^At)}, where the plus (minus) sign is taken 
when the fiavor of one B meson is opposite to (the same as) the other. 

At KEKB, the T(45') resonance is produced with a boost of (3^ — 0.425 nearly along the 
z axis defined as anti-parallel to the positron beam direction, and At can be determined 
as At ~ Az/{[3'^c), where Az is the z distance between the fcp and /tag decay vertices, 
Az = zcp — Zta,g. The average value of Az is approximately 200 /xm. 

The Belle detector [9] is a large-solid-angle spectrometer that includes a silicon vertex 
detector (SVD), a central drift chamber (CDC), an array of aerogel threshold Cerenkov coun- 
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tcrs (ACC), time-of-flight (TOF) scintillation counters, and an electromagnetic calorimeter 
(ECL) comprised of Csl(Tl) crystals, all located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron fiux-return located outside of the coil is instru- 
mented to detect mesons and to identify muons (KLM). 

II. EVENT SELECTION AND RECONSTRUCTION 

A. Reconstruction of — >^ charmonium K^*^^ decays 

We reconstruct decays to the following CP eigenstates [11]: J/ipKg, ip{2S)Kg, Xci^s^ 
TjcK^ for if = -1 and J/^jKI for = +1. We also use S° Jj^K*^ decays with the 
subsequent decay K""^ — > /T^tt^. Here the final state is a mixture of even and odd CP, 
depending on the relative orbital angular momentum of the J/-?/' and fT*". We find that the 
final state is primarily = +1; the ^/ = —1 fraction is 0.19 ± 0.02(stat) ± 0.03(syst) [12]. 

The reconstruction and selection criteria for all jcp channels used in the measurement are 
described in detail elsewhere [3]. Jj^) and iii^S) mesons are reconstructed via their decays 
to t^t~ [i = fi,e). The ipi^S) is also reconstructed via J/V'tt+tt", and the Xci via Jji)^- 
The ?7c is detected in the K'^K'n^ , K^K~ti^ , and pp modes. For the J/ipKg mode, we use 
Kg — >• 7r+7r~ and 7r°7r° decays; for other modes we only use Kg — >• tt~^ti~. For reconstructed 
B — > fcp candidates other than J/i/jK^, we identify B decays using the energy difference 

AE = - E^'^^^ and the beam-energy constrained mass Mbc = ^/(E^J^^^Jpf^ , 
where -Ebcam ^^'^ beam energy in the center-of-mass system (cms) of the T(4S') resonance, 
and E"^^^ and are the cms energy and momentum of the reconstructed B candidate, 
respectively. 

Candidate B^ J/i/jK^^ decays are selected by requiring ECL and/or KLM hit patterns 
that are consistent with the presence of a shower induced by a Kl meson. The centroid of 
the shower is required to be within a 45° cone centered on the K^ direction inferred from 
two-body decay kinematics and the measured four-momentum of the J/ip. 

B. Reconstruction of flavor-eigenstate samples 

1. B° ^ D*-e+iy 

We use the decay chain S° D*"£+z/, D*~ D^n~ , where D° K^n~, K^ti'ti^ or 
K^-K^-K^-K^ . We require associated SVD hits and radial impact parameters dr < 0.2 cm 
for all tracks. Track momenta in the laboratory frame for — > i^+7r~7r^7r~ decays are 
required to be larger than 0.2 GeV/c, while no additional requirements are applied for the 
other modes. Charged kaons arc identified by combining information from the TOF, ACC 
and dE/dx measurements in the CDC. Photon candidates are defined as isolated ECL clus- 
ters of more than 20 MeV that are not matched to any charged track, tt" candidates are 
reconstructed from pairs of photon candidates with invariant masses between 124 and 146 
MeV/c^. A mass-constrained fit is performed to improve the 7r° momentum resolution. A 
minimum 7r° momentum of 0.2 GeV/c is required. For K^tt~ and K^-k^-k^ti' can- 

didates, we use daughter combinations that have an invariant mass within 0.013 GeV/c^ of 
rriDO] for K^ir^n^ we expand the mass window to —0.037 GeV/c^ and +0.023 GeV/c^. 

For D*" — > D'^tt" decays, we combine D° candidates with a low-momentum tt" (slow pion) 
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that is reconstructed using a vertex constraint and require the mass difference between the 
D*^ and candidates, Mdig, to be within 1 MeV/c^ of the nominal value. We reject D*~ 
candidates with cms momentum greater than 2.6 GeV/c, which is beyond the kinematic 
limit for B meson decays. 

For the associated lepton, we use electrons or muons with a charge opposite to that 
of the D*~ candidate. Electron identification is based on a combination of CDC dE/dx 
information, the ACC response, and the energy deposition of the associated ECL shower. 
Muons are identified by comparing information from the KLM to extrapolated charged 
particle trajectories. We require 1.4 GeV/c < pf^ < 2.4 GeV/c, where pf^ is the cms 
momentum of the lepton. The cms angle of the lepton with respect to the direction of 
the D*~ candidate is also required to be greater than 90 degrees. For — > D*^i^u 
decays, the energies and momenta of the B meson and the D*i system in the cms satisfy 
M2 = {Ef' - Ef^lf - Ipi'^'^P - IpbTP + 2|Pb°'1 \Pdh\ cos where is the neutrino 

mass and 6b,d*£ is the angle between p^™*^ and p^T^. We calculate cos Ob,d*£ setting = 0. 
The signal region is defined as \ cos6B,D*e\ < l-l- We also require the candidate D*~i^i' 
decays to be outside the signal region when we artificially reverse the lepton momentum 
vector; this ensures that entries in the signal region after the reversal, which are used for 
the background estimation, do not contain D*~£'^i/ events. 

2. Hadronic modes 

B^ and B^ mesons are fully reconstructed in the following decay modes: D~7r~^, 
D*-7i+, D*-p+, J/ipK*^, B+ Wn+, and J/ipK+. We also use 5° ^ J/ipK'^s ^ecays, 
which are described in the previous section, with no flavor assignment. D*~ and D° can- 
didates are reconstructed in the same decay modes that are used for the D*~i'^i/ mode. 
Charged D candidates are reconstructed in the D~ K^ti~tt~ channel. For D and D*~ 
candidates, we apply mode-dependent requirements on the reconstructed D mass (ranging 
from ±15 MeV/c^ to ±50 MeV/c^) and the Mdig (ranging from ±3 MeV/c^ to ±12 MeV/c^), 
in a similar way as for the D*~i'^i' mode. Candidate K*^ — > X+tt" decays are required to 
have an invariant mass within 75 MeV/c^ of the nominal K*^ mass, p'^ candidates are 
selected as Tr+vr" pairs having invariant masses within 150 MeV/c^ of the nominal p"*" mass. 

To reduce background from e+e" qq {q = u, d, s or c) continuum events, a selection 
based on the ratio of the second to zeroth Fox- Wolfram moments [13] and the angle between 
the thrust axes of the reconstructed and associated B mesons is applied mode by mode. 

C. Flavor tagging 

For neutral B to fcp and /av decays, charged leptons, pions, kaons, and A baryons 
that are not associated with the reconstructed decay are used to identify the 6-fiavor of the 
accompanying B meson. The tagging algorithm is described in detail elsewhere [14]. We 
use two parameters, q and r, to represent the tagging information. The first, q, has the 
discrete value +1 (—1) when the tag-side B meson is likely to be a B^ (-8°)- The parameter 
r corresponds to an event-by-event flavor-tagging dilution that ranges from r = for no 
flavor discrimination to r = 1 for an unambiguous flavor assignment. It is determined from 
a large number of events generated by Monte Carlo (MC) simulation, and is used only to 
sort data into six intervals of r, according to estimated flavor purity. We determine directly 
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from data the average wrong-tag probabilities, wi = (w^' + tVi )/2 (Z = 1, 6), and differences 
between and decays, Awi = Wi~ — wf, where w^^'"^ is the wrong-tag probabihty for 
the B^[B^) decay in each r intervaL 

D. Vertex reconstruction 

The vertex position for the fcp decay is reconstructed using leptons from J/ip decays or 
charged hadrons from rjc decays. Each vertex position is required to be consistent with the 
interaction-region profile (IP), determined run-by-run, smeared in the r-0 plane to account 
for the B meson decay length. With the IP constraint, we are able to determine a vertex 
even with a single track; the fraction of single-track vertices is about 10% for zcp- 

The vertex position for the /av decay that includes a D meson in its decay products is 
reconstructed using the D meson trajectory, a track other than the slow vr" candidate from 
D*^ decay, and the IP constraint. For the J/ipK^ and B^ J /iI)K*^{K^ti^) decays, 

we use leptons from J/ if) decays and the IP constraint, in the same way as for the fcp decay. 

The vertex position for the /tag is obtained with the IP constraint and with well recon- 
structed tracks that are not assigned to fcp or /av. The algorithm is described in detail 
elsewhere [15]. The fraction of single-track vertices is about 22% for 2;tag. 

We only use events with vertices that satisfy |At| < 70 ps and ^ < 100, where ^ is the 

of the vertex fit calculated only in the z direction [15]. The overall vertex reconstruction 
efficiency is 87.1 ± 0.7% for B^ J/i^K^s candidates. 

The proper-time interval resolution function Rsig{At) is formed by convolving four com- 
ponents: the detector resolutions for zcp and Ztag, the shift in the 2;tag vertex position due 
to secondary tracks originating from charmed particle decays, and the kinematic approxi- 
mation that the B mesons are at rest in the cms [15]. A small component of broad outhers 
in the Az distribution, caused by mis-reconstruction, is represented by a Gaussian function. 

E. Signal yields 

After flavor tagging and vertexing, we find 5417 fcp candidates in total in the signal 
region; these are used for the sin 20i determination. Table I lists the numbers of candidates, 
A^ev; a-nd the estimated signal purity for each fcp mode. Figure 1 shows the Mbc distribution 
after applying mode-dependent requirements on AE for all B^ candidates except for B^ 
J/ipK^. There are 3085 entries in total in the signal region defined as 5.27 GeV/c^ < Mbc < 
5.29 GeV/c^. Figure 2 shows the pg"*^ distribution for B^ — > J/ipK^ candidates. We find 
2332 entries in the 0.20 GeV/c < p''^ < 0.45 GeV/c signal region. 

Table II lists A'cv and the purity for each /av mode after the vertexing. The total number 
of /flv candidates is 177368 with a purity of 81%. Figure 3 shows the cos9b,d*i distribution 
for the D*~i^i' candidates. Figure 4 shows the Mbc distributions for B^ and B^ decays to 
/flv states. 

According to a MC simulation study, there is a small fraction of background (less than 
1% for fcp and 3% for /gy candidates) from other B decays peaking in the Mbc signal region. 
The effect of the peaking background is treated as a systematic error. 
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TABLE I: Numbers of reconstructed B — > fcp candidates after flavor tag 


ging and vertex recon- 


srrucxion, iVev) ana tne esximateu 


signal purity, p. 








Mode 






A^ev 


P 




-1 




1997 


0.976 ± 0.001 




-1 




288 


0.82 ±0.02 




-1 




145 


0.93 ±0.01 


i){2S){J/i)TT+TT-)Kl 


-1 




163 


0.88 ±0.01 




-1 




101 


0.92 ±0.01 




—1 




123 


0.72 ±0.03 




-1 




74 


0.70 ±0.04 




-1 




20 


0.91 ±0.02 


All with = -1 


-1 




2911 


0.933 ± 0.002 




+1(81%) 




174 


0.93 ±0.01 


J/i^Kl 


+1 




2332 


0.63 ±0.03 


TABLE IL Numbers of reconstructed B ^ /gy candidates 


after vertex reconstruction, Aev, and 


the estimated signal purity, p. J/ipKg candidates 


are used with no flavor assignment. 


Mode 




iVev 




P 


D*-e+iy 




84823 




0.781 






11921 




0.888 






11156 




0.899 






8767 




0.763 


J/iljK*^{K+n-) 




3681 




0.954 


J/i>K^{Tr+Tr-) 




2001 




0.976 


5° total 




122349 




0.809 






46248 




0.783 


J/tpK+ 




8771 




0.966 


B+ total 




55019 




0.812 


B^ + B+ total 




177368 




0.810 



III. FIT RESULTS WITH FLAVOR-EIGENSTATE SAMPLES 

We perform a multi-parameter fit to fiavor-eigenstate samples to obtain the B meson 
lifetimes, the B'^-B'^ mixing parameter Am^, wrong-tag fractions, and parameters for the 
resolution function simultaneously. We use 5° — > D*~t^v, J/iIjK*^(K~^7t~), D*~7r+, D~7r~^, 
D*~p'^, and J/iljKg{i^i~) (for t^o and resolution parameters only) for decays, and 
B~^ — > D'^TT^ and J/ipK~^ for B~^ decays. The fit uses 32 parameters; 12 for wrong-tag 
fractions, 14 for the resolution function, 3 for the S+ background in B^ decays, and 3 
physics parameters Am^, r^o and tb+- We also obtain the lifetime ratio, r^-^ = tb+/tbo, 
by repeating the fit in which tb+ is replaced with Tt-^Tbo. Two of the 14 parameters for the 
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FIG. 1: Beam-energy constrained mass distribution within the AE signal region for all fcp modes 
other than J/ipK^. The solid curve shows the fit to signal plus background distributions, and the 
dashed curve shows the background contribution. 



resolution function are newly added to the resolution function described in [15] to improve 
the description of the effect of charmed particle decays on the /tag vertex. 

The probability density function (PDF) expected for the signal distribution for B° decays 
to /flv is given by 

V^r\At, q, wi, Awi) = g^^^^ {l - qAwi + [-](1 - 2wi) cos(Am,Ai)}, (2) 

where OF (SF) denotes _B°_B° {B^B'^ or B^B^), i.e. a state with the opposite (same) flavor. 
The signal PDF for 5+ decays is given by 



^, e-\^t\/rB+ 



(3) 



The signal PDFs are convolved with the i?sig(At) to account for the detector resolution. 

The background PDF for the hadronic modes is modeled as a sum of exponential and 
prompt components, 

-|At|/Tbkg 

nkg(At) = (1 - fs)— + f55{At), (4) 

ZTbkg 

(1— Z^) is the fraction of the exponential component with the effective lifetime Tbkg, and 5{At) 
is the Dirac delta function. It is convolved with a sum of two Gaussians, which is used as the 
background resolution function. The parameters for the background PDF are determined 
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FIG. 2: p^j^^ distribution for B'^ J/t/jK^ candidates with the results of the fit. The dashed Unes 
indicate the signal region. 
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FIG. 3: cos9B,D*e distribution for the D*~t^u candidates. The circles with errors show the data. 
The solid line is the fit result. The total background and the D**ii' component are shown by the 
dashed line and the hatched area, respectively. The inset shows the same figure with a logarithmic 
vertical scale. 
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FIG. 4: Beam-energy constrained mass distributions in the /S.E signal region for flavor-eigenstate 

candidates in neutral B meson decays (left) and in charged B meson decays (right). The solid curves 
show the fits to signal plus background distributions and the dashed curves show the background 
contributions. 



using the A£'-Mbc sideband region for each decay mode. For D^n'^ decays, using 

events outside the signal region, the value for fs is determined to be 0.49 ±0.01 (0.45 ±0.03) 
for events with multi-track (single-track) vertices and the effective lifetime Tbkg is found to 
be 0.93 ± 0.03 ps. The parameters for other B f^^ decays that include a D meson as 
a decay product arc similar to those for the B^ — >■ D^tt^ decay. A similar procedure for 
B+ J/i/jK+ decays yields fs = 0.86 ± 0.05 (0.75 ± 0.08) for events with multi-track 
(single-track) vertices and Tbkg = 1-55 ± 0.22 ps. The parameters for B^ — > J/ijjKg are 
similar to those for B^ J/ipK^. The value for fs in B^ J /iI)K*^{K^tt~) decays is 
found to be small; fs is 0.07 ± 0.05 for events with multi-track vertices and is fixed at zero 
for events with single-track vertices. The effective hfetime Tbkg is 1.50 ± 0.05 ps. 

The background for the D*~l^v decay is divided into four components: B — > £)**£z/ events 
(8.7 ±0.3%); fake D* mesons (8.0 ±0.1%); random combination of D* mesons with leptons 
with no angular correlation called "uncorrelated background" (2.5±0.1%); continuum events 
(2.7 ± 0.2%). Here D** consists of charmed mesons heavier than the D* meson and non- 
resonant D*7r components. The PDF for the B — > D**iu background is given by a sum of 
B° and S"*" components, 

^"^^(A*) = (1 - MPZ""' + fB.V°/:'-\ (5) 

where /b+ is the B^ fraction in the B — > D**ii' background. The is given by 

pOF(At) = (1 _ w'^+)P^^+{At) and Vf+{At) = w'j,+P^^+{At), where is the wrong tag 
fraction determined from the B~^ D^tt"*" sample and is given by 

<(At) - (1 - fr'J^^ + (6) 



B+ 
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Here /,-' and are the fraction and the effective hfetime for events in which an additional 
TT'^ from the D** decays contaminates the /tag vertex reconstruction. The parameters /^+, 
/t-' and t'q+ arc determined in the final fit. To determine these parameters precisely, 
events in —10 < cos9B,D*i < —1-1, where the D**ii' background events are dominant, 
are also included in the fit. The fit yields /b+ = 0.51 ± 0.04, = 0.56 ± 0.10 and 

= 0.74 ±0.14 ps. 

For continuum and uncorrelated backgrounds, the same functional form as that of the 
hadronic background PDF is used. The parameters for continuum are determined from off- 
resonance data to be fs = 0.55 ±0.09 (0.58 ±0.11) for events with multi-track (single-track) 
vertices and rbkg = 0.80 ± 0.08 ps. For the uncorrelated background, a fit to the sample 
outside the signal region yields = 0.15 ± 0.08 and Tbkg = 1-23 ± 0.06 ps. 

The PDF of the fake D* background is given by Eq.(4) with a mixing component added 
to account for oscillation in the background. A fit to events in the M^is sideband yields 
fs = 0.13 ± 0.02 (0.04 ± 0.03) for events with multi-track (single-track) vertices, Tbkg = 
1.49±0.03 ps and Am^kg — 0.54±0.05 ps~^, where Ambfcp is the effective mixing parameter. 
The fraction of the mixing component and wrong tag fractions are determined for each of 
the six intervals of the flavor tag quality r. The wrong tag fractions range from 0.50 ± 0.01 
for the lowest r region to 0.20 ± 0.01 for the highest r region. The fraction of the mixing 
component for the lowest r region is fixed at 0. Values for the other r intervals range from 
0.39 ±0.12 to 0.84 ±0.09. 

We test the fit method and parameterization with a large number of MC events, and 
obtain results consistent with the input values. The wrong-tag fractions obtained with the 
MC events are also found to be correct. 

The unbinned maximum-likelihood fit to data yields 



The results are consistent with our previous measurements [6-8] and supersede them. 
Figure 5 shows the flavor asymmetry, A{At) = [NoF(At) - NsF(At)]/[NoF(At) + NsF(At)], 
where Nof(sf) denotes the number of OF (SF) events. The results of the lifetime measure- 
ments for neutral and charged B meson decays are shown in Fig. 6. 

Systematic uncertainties are listed in Table III. The method to determine the systematic 
errors due to the vertex reconstruction follows the same procedure as for the sin 201 measure- 
ment, which will be explained later. We estimate the contribution due to uncertainties in 
the resolution function by comparison of different parameterizations, as well as by changing 
parameters that are derived from MC to model the effect of non-primary tracks [15]. A pos- 
sible bias in the event reconstruction and fitting procedure is checked with a large number 
of MC events. We find no bias and take the statistical error in MC as a systematic error. 
Several D** components are used in this analysis to model the cos 9b,dh shape for the D**iu 
background. To estimate the systematic errors due to uncertainties of the fractions of the 
D** components, we set the fraction of each component to unity (with all other components 
set to zero) and repeat the analysis; for each measurement, we take the largest variation 
on the result as the systematic error. Systematic errors that arise from uncertainties in 
other background fractions and from the background At shape are obtained by varying each 



rB+ 

Tb+Z^bo 

Arud 



[1.534 ± 0.008(stat) ± O.OlO(syst)] ps, 
[1.635 ± O.Oll(stat) ± O.Oll(syst)] ps, 
1.066 ± 0.008(stat) ± 0.008(syst), 



[0.511 ± 0.005(stat) ± 0.006(syst)] ps 
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FIG. 5: Time-dependent flavor asymmetry for flavor-eigenstate decays. The curve is the result of 
the unbinned maximum-likelihood fit. 
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FIG. 6: At distributions of neutral B meson pairs (left) and charged B meson pairs (right). 
The solid lines represent the results of the unbinned maximum- likelihood fit. The dashed lines 
correspond to the outlier components. 
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TABLE III: Summary of the systematic errors on the measurement of r^o , Tb+ , / t^o and Am^. 



Source 


TbO 




tb+/tbo 


Airid 


Vertex reconstruction 


0.005 


0.007 


0.003 


0.003 


Resolution function 


0.004 


0.005 


0.004 


0.001 


Possible fit bias 


0.003 


0.004 


0.003 


0.002 


D**iu background 


0.004 


0.003 


0.002 


0.004 


Other background fraction 


0.001 


0.006 


0.003 


0.001 


Background At shape 


0.006 


0.003 


0.005 


0.002 


Total 


0.010 


0.011 


0.008 


0.006 




(ps) 


(ps) 




(ps-^ 



TABLE IV: Event fractions ei, wrong-tag fractions wi, wrong-tag fraction differences Awi, and 
average effective tagging efficiencies egg = e/(l — 2wi)'^ for each r interval. The errors include both 
statistical and systematic uncertainties. 



I 


r interval 


ei 


Wl 


Awi 




1 


0.000 


- 0.250 


0.398 


0.464 ± 0.006 


-0.011 ±0.006 


0.002 ± 0.001 


2 


0.250 


- 0.500 


0.146 


0.331 ± 0.008 


+0.004 ± 0.010 


0.017 ±0.002 


3 


0.500 


- 0.625 


0.104 


0.231 ± 0.009 


-0.011 ±0.010 


0.030 ± 0.002 


4 


0.625 


- 0.750 


0.122 


0.163 ± 0.008 


-0.007 ± 0.009 


0.055 ± 0.003 


5 


0.750 


- 0.875 


0.094 


0.109 ± 0.007 


+0.016 ± 0.009 


0.057 ± 0.002 


6 


0.875 


- 1.000 


0.136 


0.020 ± 0.005 


+0.003 ± 0.006 


0.126 ±0.003 



parameter individually, repeating the fit procedure, and adding each contribution in quadra- 
ture. In the nominal fit, we do not include a mixing component in the background PDF for 
the hadronic decays. We repeat the fit with a background PDF including a mixing term. 
Uncertainties in the overall z scale of the detector arising from the measurement error and 
thermal expansion during the operation are found to be negligible. 

The same fit also yields wrong-tag fractions that are summarized in Table IV. The total 
effective tagging efficiency is determined to be e^s = Z^Li ^((1 ~ 2w;)^ = 0.287 ± 0.005, 
where ei is the event fraction for each r interval determined from the J/ ipKg simulation and 
is listed in Table IV. The error includes both statistical and systematic uncertainties. 

We find that the average At resolution is ~ 1.43 ps (rms). The width of the outher 
component is determined to be (39 ± 2) ps; the fractions of the outlier components are 
(2.1 ±0.6) X 10~^ for events with both vertices reconstructed with more than one track, and 
(3.1 ±0.1) X 10~^ for events with at least one single-track vertex. 

IV. RESULTS OF CP ASYMMETRY MEASUREMENTS 

Figure 7 shows the observed At distributions for the q^j = +1 and q^f = —1 event 
samples (top), the asymmetry between two samples with < r < 0.5 (middle) and with 
0.5 < r < 1.0 (bottom). The asymmetry in the region 0.5 < r < 1.0, where wrong-tag 
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FIG. 7: At distributions for the events with q^f = —1 (open points) and q^f = +1 (sohd points) 
with ah modes combined (top), asymmetry between g^j = —1 and qif = +1 samples with < 
r < 0.5 (middle), and with 0.5 < r < 1 (bottom). The results of the global unbinned maximum- 
likelihood fit (sin 201 = 0.728) are also shown. 



fractions are small as shown in Table IV, clearly demonstrates large CP violation. 

We determine sin 20i from an unbinned maximum-likelihood fit to the observed At dis- 
tributions. The PDF for the signal distribution is given by 



VsigiAt, q,wi,Awi,^f) 



-,-\At\/TgO r 

Atbo 



qAwi — qC,f{l — 2wi)sm2(l)i sm{AmdAt) 



where we fix the lifetime r^o and mass difference Am^^ at their world average values [16]. 
Each PDF is convolved with the appropriate i?sig(At) to determine the likelihood value for 
each event as a function of sin20i: 



ol 



+ 00 



/sigPsig(At')i?sig(At - At') 
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+ (1 - /sig)nkg(AOi?bkg(Ai - Ai') dAt' + UP,,{At), 



(12) 



where /sig is the signal fraction calculated as a function of for J/ipK^ and of AE and 
Mbc for other modes. Pbkg(At) is the PDF for combinatorial background events, which is 
modeled as a sum of exponential and prompt components. It is convolved with a sum of 
two Gaussians, i?bkg, which is used as the background resolution function. We assume no 
asymmetry in the background At distribution. 

For S° —>■ J/ipKi^ and J/ipK*^ decays, in addition to the combinatorial background, 
background events from other B decays and their CP asymmetries are considered. The 
background in the J/ipKl mode is dominated by the following B J/i/jX decays: J/ijjKg 
having = -1 (10 ± 2%); ij{2S)Kl, xaKl and J/^7r° having if = +1 (4 ± 1%); 
J/^K*\K*^ KIti^) (20 ± 2%), which is a mixture of = -1 (81%) and if = +1 
(19%); other non-CP modes (66 ± 1%). The fraction of each component is obtained from 
a MC simulation study. For J / ip K*^ (K^tt^) decays, we include in our PDF contaminations 
from other B — > J/ipK* decays (7.1 ± 0.3%) and non-resonant B^ J/ipK^n^ decays 
(6.3 ± 0.5%) in the Mbc peak. The background fractions are obtained from MC and the K* 
mass sideband [12]. We use the signal PDF with no CP asymmetry for these components. 

To account for a small number of events that give large At in both the signal and back- 
ground, we introduce the PDF of the outlier component, Poh ^-iid its fraction /oi. The only 
free parameter in the final fit is sin20i, which is determined by maximizing the likelihood 
function L = JJiPi, where the product is over all events. We obtain 



The result is consistent with the value in our previous publication [4] and supersedes it with 
a reduced error. 

The signal PDF for a neutral B meson decaying into a CP eigenstate [Eq. (11)] can be 
expressed in a more general form as 



where S = 21m(A)/(|Ap + 1), ^ = (|Ap — l)/(|Ap + 1), and A is a complex parameter 
that depends on both B^-B^ mixing and on the amplitudes for B^ and B^ decay to a CP 
eigenstate. The presence of the cosine term (|A| ^ 1) would indicate direct CP violation; the 
value for sin 20i reported above is determined with the assumption | A| = 1, as | A| is expected 
to be very close to one in the SM. In order to test this assumption, we also performed a fit 
using the expression above with acp = — ^/Im(A)/|A| and |A| as free parameters, keeping 
everything else the same. We obtain 



and acp — 0.728 ± 0.056(stat). This result is consistent with the assumption of no direct 
CP violation used in our analysis and the acp term is in good agreement with the sin20i 
value obtained with the one-parameter fit. 

Table V lists the systematic errors on sin20i and |A|. The total systematic error is 
obtained by adding each of them in quadrature. The largest contribution for sin 201 comes 



sin 201 = 0.728 ± 0.056(stat) ± 0.023(syst). 



(13) 




(14) 



A| = 1.007 ± 0.041(stat) ± 0.033(syst), 



(15) 
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TABLE V: Summary of the systematic errors on sin2(/>i and |A|. 



Source 


sin 2<f)i 




Vertex reconstruction 


0.013 


0.012 


Th I 51'^7"i~»T* CO"! Tier 


007 


008 


RpRolution fiinption 


0.008 


0.004 


Fossibie lit bias 


O.OOo 


O.OUd 


Background fraction (J/i/jK^) 


0.011 


0.003 


Background fraction (except for J/tpK^) 


0.007 


0.007 


Physics {tbo, Anid, J/ipK*^) 


0.003 


0.001 


Background At shape 


0.002 


0.001 


Tag-side interference 


0.002 


0.028 


Total 


0.023 


0.033 



from vertex reconstruction. The systematic error due to the IP constraint in the vertex 
reconstruction is estimated by varying (±10 /im) the smearing used to account for the B 
flight length. The track selection criteria are also varied to search for possible systematic 
biases. The effect of the vertex quality cut is estimated by varying the cut to ^ < 50 and 
^ < 200. We vary the |At| range by ±30 ps to estimate the systematic uncertainty due 
to the I At I fit range. Small biases in the Az measurement are observed in e~^e~ — > 
and other control samples. Systematic errors are estimated by applying special correction 
functions to account for the observed biases, repeating the fit, and comparing the obtained 
values with the nominal results. Systematic errors due to imperfect SVD alignment are 
determined from MC samples that have artificial mis-alignment effects to reproduce impact- 
parameter resolutions observed in data. In these studies, whenever required, we repeat the 
fit to the /flv samples, update resolution function parameters and wrong tag fractions, and 
perform the fit to CP-cigcnstatc event samples using the updated parameters so that the 
uncertainties in question arc treated in a consistent way. 

Systematic errors due to uncertainties in the wrong tag fractions given in Table IV are 
studied by varying the wrong tag fraction individually for each r region. Possible differences 
of the tagging performance between fcp and /av events are estimated using MC events. 

Systematic errors due to uncertainties in the resolution function are estimated by vary- 
ing each resolution parameter obtained from data (MC) by ±1(7 (±2(7), repeating the fit 
and adding each variation in quadrature. We also divide the entire data set into two and 
prepare two sets of resolution parameters to consider a possible difference in the detector 
performance. We repeat the fit with these resolution parameters and assign the difference 
from the nominal result as a systematic error. We also include other sources examined for 
the fit to the fiavor-eigenstate samples, which are explained in the previous section. 

A possible fit bias is examined by a fit to a large number of MC events. We find no bias 
and take the statistical error from the MC as a systematic error. 

Systematic errors from uncertainties in the background fractions and in the background 
At shape are estimated by varying each background parameter obtained from data (MC) by 
±lcr (±2cr). The systematic error due to CP content in the J/ipK^ backgrounds is checked 
by varying the parameters obtained from the MC by ±2(7. 

The small peaking background in the M^^ signal region of fcp modes other than J/ipK*^ 
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TABLE VI: Numbers of candidate 


events, 


iVev, and values of sin2^i, 


|A| for various subsamples 


(statistical errors only). 








Sample 


iVev 


sin 2^1 


1^1 


J/i)Kl{TT+TT-) 


1997 


0.67 ± 0.08 


0.98 ± 0.06 




288 


0.72 ± 0.20 


1.18 ±0.27 




308 


0.89 ± 0.20 


0.94 ±0.14 




101 


1.54 ± 0.49 


0.76 ± 0.22 




217 


1.32 ±0.28 


1.10 ±0.30 


All with = -1 


2911 


0.73 ± 0.06 


0.99 ±0.05 




2332 


0.77 ±0.13 


1.04 ±0.08 




174 


0.10 ±0.45 


1.11 ±0.33 


/tag = B^ {q = +1) 


2717 


0.72 ± 0.09 


0.89 ± 0.09 


/tag = B^ iq = -1) 


2700 


0.74 ± 0.08 


1.17 ±0.11 


< r < 0.5 


2985 


0.95 ± 0.26 


1.18 ±0.22 




1 OOA 


n (58 1 nil 


1 1 1 -1- n no 
i.ii zt u.uy 


0.75 < r < 1 


1208 


0.73 ± 0.07 


0.95 ± 0.05 


Data set I (78 fb"!) 


3013 


0.72 ± 0.07 


0.95 ± 0.05 


Data set II (62 fb"^) 


2404 


0.74 ± 0.09 


1.09 ± 0.07 


All 


5417 


0.728 ± 0.056 


1.007 ±0.041 



is neglected in the nominal analysis. The effect of the fractions and their CP asymmetries 
is studied with MC simulation and is included in systematic errors. 

Each physics parameter (tbo, Am^, J/tpK*^ polarization) is also varied by its error; for 
Am^, we also use our result (Am^ = 0.511 ps~^), repeat the fit and take the larger change 
as the systematic error. 

Finally, wc investigate the effects of interference between CKM-favored and CKM- 
suppressed B ^ D transitions in the /tag final state [17]. A small correction to the PDF for 
the signal distribution arises from the interference. We estimate the amount of correction 
using the D*~l^v sample. We then generate MC pseudo-experiments and make an 

ensemble test to obtain systematic biases in sin20i and |A|. We find that the effect on 
sin 201 is negligibly small, while a possible shift in |A| becomes the largest contribution to 
the systematic error. 

Several checks on the measurement are performed. Table VI lists the results obtained by 
applying the same analysis to various subsamples. All values are statistically consistent with 
each other. Figure 8 shows the raw asymmetries and the fit results for (cc)A'^ (top) and 
J/ipKl (bottom). A fit to the non-CP eigenstate modes B^ D* l^v and J / ipK*^[K^Ti ), 
where no asymmetry is expected, yields "sin20i"= 0.012 ± 0.013(stat). 

V. SUMMARY 

Using 152 X 10^ BB pairs collected at the T{AS) resonance with the Belle detector at the 
KEKB asymmetric-energy e'^e" collider, we have measured the CP- violation parameters 
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FIG. 8: Raw asymmetries for {cc)Kg (^j = —1) (top) and J/ifjK^ i^f — +1) (bottom). The cm'ves 
are the results of the global unbinned maximum-likelihood fit. 



sin2(/)i and |A|, -B meson lifetimes and their ratio, and the mixing parameter Arrid. These 
are basic parameters of the standard model. The results are summarized as follows: 

sin20i = 0.728 ± 0.056(stat) ±0.023(syst), 

|A| = 1.007 ± 0.041(stat) ±0.033(syst), 

tbo = [1.534 ± 0.008(stat) ±0.010(syst)] ps, 

tb+ = [1.635 ± O.Oll(stat) ±0.011(syst)] ps, 

tb+Ztbo = 1.066 ± 0.008(stat) ±0.008(syst), 

Arrid = [0.511 ± 0.005(stat) ± 0.006(syst)] ps"\ 

All results are significant improvements in precision from the previous measurements, and 
are in agreement with the standard model expectations. The significance of the observed 
deviation from unity in the lifetime ratio exceeds five standard deviations for the first time 
by a single measurement. 
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